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Baryonic matter is studied in the Skyrme model by taking into account the roles of tt, 
p, and u) mesons through the hidden local symmetry up to Oip"^) terms including the 
homogeneous Wess-Zumino (hWZ) terms. Using the master formulas for the low energy 
constants derived from holographic QCD models the skyrmion matter properties can be 
^ ^ quantitatively calculated with the input values of the pion decay constant /,r and the 

On ■ vector meson mass m . We find that the hWZ terms are responsible for the repulsive 

QQ ' interactions of the ui meson. In addition, the self-consistently included 0{p'^) terms with 

^v^ ' the hWZ terms is found to increase the half skyrmion phase transition point above the 

f**^ ^ normal nucleon density. 

^^ , Keywords: Baryonic matter; Skyrme model; Hidden local symmetry 

o 

1. Introduction 

. ^ ^ Understanding the QCD phase structure is an extremely important problem in 

j^ \ particle and nucleon physics, but is highly nontrivial because of the nonperturbative 

H ' nature of strong interactions. Therefore, any theoretical investigation of the QCD 

- - - phase structure requires effective models for strong interactions at low energies. In 

this talk, we discuss the QCD phase structure for zero temperature baryonic matter 
in the Skyrme model based on the hidden local symmetry (HLS) up to 0{p'^) terms 
including the homogeneous Wess-Zumino (hWZ) terms.'i''^ To uniformly describe 
the baryon, baryonic matter and also the change of hadron properties in baryonic 
matter background, we apply the Skyrme's conjecture^ and describe baryons as 
topological solitons in a mesonic theory. We also put the skyrmions into the face- 
centered-cubic (FCC) crystaP for studying baryonic matter. 

Although the original Skyrme model is quite simple as it is just the nonlinear 
sigma model stabilized by the Skyrme term, it can reasonably describe the prop- 
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erties of a single nucleon in vacuum in the sense of large Nc expansion.® Since the 
original Skyrme model contains only pion degree of freedom, it is natural to improve 
the model predictions by including more dynamical degrees of freedom. When we 
increase the energy scale, the first resonance entering is the ground state p and w 
vector mesons. In the present work, we adopt the Hidden Local Symmetry (HLS) 
approach which is gauge equivalent to the nonlinear sigma model for including vec- 
tor mesons.!^ 

Once extra degrees of freedom are added to an effective theory, the dynamics 
generally becomes more complicated by introducing more independent low energy 
constants (LECs) which are to be determined by empirical data. However, it is 
not possible to determine all the LECs phenomenologically without ambiguity as 
the number of LECs increases, and this hinders the systematic investigation of the 
baryon and baryonic matter properties. In Ref.Elwe derived general master formulas 
for the LECs of HLS from a kind of holographic QCD (hQCD) models. This allows 
us to investigate the skyrmion and hadronic matter properties with only two input 
values of /,r and rup. 

In Refs. I1I2I we applied the HLS up to 0{p'^) terms including the hWZ terms 
to investigate the single-skyrmion properties. It was found that the hWZ terms are 
responsible for the role of the u meson and that the w and p mesons have very 
different roles indicated by the attractive nature of the p meson and the repulsive 
interaction from the cj meson. As the next step for studying the role of vector 
mesons, we apply the above model to baryonic matted to study the skyrmion 
matter properties and the half skyrmion phase transition point. 

This paper is organized as follows. In Sec. [2] we investigate the single skyrmion 
properties from the HLS Lagrangian after discussing the implications of the infinite 
towers of vector mesons in the skyrmion properties. We then study the skyrmion 
matter properties from HLS in Sec.[3l where we also compare our results with other 
model calculations. Section |4] contains a summary and perspective of our work. 

2. Skyrmion properties from the hidden local symmetry 

We start with the study of the single skyrmion properties based on the HLS up to 
O(p^) terms including the hWZ terms. Therefore, this model includes the it, p and 
ljj mesons. The LECs of the HLS Lagrangian can be determined from the master 
formulas derived from hQCD models which allows us to control the LECs so that 
the skyrmion properties can be calculated in a self-consistent way. 

The role of the infinite tower structure in five-dimension hQCD models in the 
skyrmion properties have been addressed in the large 't Hooft coupling A limit.'^It 
was found that, when only the pion was kept in the truncated model, the skyrmion 
has an excessive energy of about 25% of the Bogomol'nyi bound, while the excess 
energy is reduced to 7% when the p meson is included. This observation indicates 
that high-lying vector mesons make the theory flow to a conformal theory. However, 
in the large A limit, the theory does not include the repulsive interactions due to the 
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w meson which is the corrections of the next order in A, where the Chern-Simons 
(CS) term brings in a U(l) degree of freedom responsible for the w meson and its 
tower. 

To include the u) meson, we consider the hWZ terms in the HLS Lagrangian, 
which allows us to study the role of the ui meson in the baryon properties. In 
addition, since the validity of the Skyrme model is supported by the large Nc limit, 
for consistency, we need to take into account the higher order terms such as the 
0{p^) terms in addition to the 0{p^) terms. Our Lagrangian is then summarized as 

-ChLS — -^(2) + -^(4) + -^anom, (1) 

where £(2), -^(4) and £anom stand for the O(p^), C(p'*) and the hWZ terms, respec- 
tively. The explicit expressions for the Lagrangian ([T]) responsible for the calculation 
of the skyrmion properties can be found in Ref. JL 

The Lagrangian ^ has seventeen LECs, namely, /^, 5, a, yi ^ yg, Z4, Z5 and 
ci ~ C3. These LECs can be determined from the general master formulgpJ^I derived 
from a class of hQCD models which is defined by the following general 5D action: 

S,^S^^' + S^^ (2) 

where the Dirac-Born-Infeld part 5*^^^ and the CS part S^^ are expressed as 

where Gym = A/(1087r^) and J^mn ~ OmAn — OnAm — «[^m,^7v] is the field 
strength of the 5D gauge field Am- Here, we use the index M — {fi,z) with 
the convention /i = 0,1,2,3 and Ki^2{z) are the metric functions constrained by 
the gauge/gravity duality. In Eq. (JH), M^ and R stand for the four-dimensional 
Minkowski space-time and z-coordinate space, respectively, and w^{A) is the CS 
5-form. Note that the DBI part is of 0{X^) while the CS term is of 0(A°). From the 
action ([2]), after integrating out all the modes except the zero and the lowest-lying 
vector modes, the LECs of HLS are expressed as 

/^ = N,Gy^M],j, I dzK^iz) [Mz)]\ II = N^G^^Ml^XM), 

\ = N.GyMi^l), y, = -2/2 = -NcGyu ((l + '/'i - ^0)' 

Vs = -Vi - -NcGym (y^i (1 + ^1)') , ^5 = 2j/8 = -2/9 = -2N,Gyu (V-^^o) 

% = - (% + Vt) > Vt - 27V,Gym {i^i (1 + V^i) (1 + ^1 - ^0)) , 

Z4 = 2N,Gyu (V-i (1 + V^i - ^0)) , ^5 - -2^cGym (V'? (1 + V-i)) , 

• /I 2 1 2 1 

V'oV'i ( 2^°"*' 6^^ ~ 2 
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where Ai is the smallest (non-zero) eigenvalue of the eigenvalue equation, 

- K-\z)d., [K2{z)d.^„{z)] = A„^„(z), (6) 

and 

/•CXD /'OO 

{A) = / dzK,{z)A{z), {{A)) = / dzA{z). (7) 



With the help of the formulas in Eq. ([5]) , sixteen of the seventeen parameters can 
be fixed from the two parameters Gym and Mkk which, in turn, can be determined 
by the pion decay constant f^^ and the vector meson mass rup. The only undeter- 
mined parameter is the HLS-parameter a or equivalcntly the coupling constant g 
as they are connected by the mass relation m'^ — ag^f^. From the point view of the 
5D holographic model, this ambiguity is attributed to the normalization of the 5D 
wave function ipiiz) that cannot be determined from the homogeneous eigenvalue 
equation ([5]). In Refs. .l;2. we have shown that the parameter g disappears from the 
expressions for the skyrmion mass if we express the LECs in terms of the normalized 
tjji{z), which leads to the independence of the skyrmion properties on the value of 
the HLS parameter a. This conclusion was numerically verified in Ref. 2, 

For numerical calculations, we use the empirical values, rUp — 775.5 MeV and 
/tt = 92.4 MeV as input values. And for the hQCD model, we take the Sakai- 
Sugimoto model,'^where Ki{z) = {l + z'^)~^/^ and K2{z) = l + z'^. Here, we consider 
three versions of the HLS model: HLS(7r, p, uj), HLS(7r, p) and HLS(7r). HLS(7r, p, w) 
is the full model including the pion, p, and lo mesons and HLS(7r, p) is the model 
without the hWZ terms. When the vector fields are integrated out, we have the 
original Skyrme model but with the Skyrme parameter determined by the hQCD 
model. This is called HLS(7r). Our numerical results are summarized in Table[T]that 
shows the skyrmion mass Mgoi, the A-A^ mass splitting Am, the winding number 
rms radius ^{r'^)w and the energy rms radius \/ {r'^)E- For a comparison, we also 
show the results of "the minimal model" which includes the uj meson in a minimal 
wajWJ in the column of HLSniin(7r, p, w). 

Table 1. Skyrmion properties calculated from the HLS models. The 
soliton mass M^^l and the A-iV mass difference Ajv/ are in units of 



MeV, while ^ {r^)^^ and ^(^ 





HLS(7r,p,aj) 


HLS(7r,p) 


HLS(7r) 


HLS 


miii(vr,p, <^) 




1184 
448 


834 
1707 


922 
1014 




1407 
259 




0.433 
0.608 


0.247 
0.371 


0.309 
0.417 




0.540 
0.725 
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Fig. 1. The energy per baryon number E jB and the expectation value (tr) of the mininiuni energy 
configuration for a given value of L. 



Close inspection of the results given in Table [T] leads to the following observa- 
tions. Namely, the results from HLS(7r) and HLS(7r,/9) explicitly show the attrac- 
tive interaction nature arising from the p meson, while comparison of the results 
of HLS(7r,p) and HLS(7r,p,w) explicitly illustrates the repulsive nature due to the 
w meson. In addition, the differences between the results from HLS(7r,p, w) and 
HLSinin(7r, p, w) come from the ©(p"*) terms and the extra terms in the hWZ terms. 



3. Skyrmion matter from the hidden local symmetry 

In this Section, we apply the HLS model of the previous Section for the study of the 
skyrmion matter properties- using the procedure explained in Ref.[5j i.e., by putting 
skyrmions into the FCC crystal. Shown in Fig.[T]is the numerical results for EjB and 
{a) as functions of the FCC size parameter L. Here, (cr) is the average of the chiral 
field V over the space where a single skyrmion occupies, which vanishes in the half- 
skyrmion phase. Here, we present the numerical results for HLS(7r,p, w), HLS(7r,p) 
and HLSi„in(7r,p, w) by the solid, dahed, and dash-dotted lines, respectively. The 
position of the normal nuclear matter density is denoted by the vertical dotted line. 

One can sec that, in HLS(7r, p, w), the half-skyrmion phase transition happens al- 
most at the normal nuclear matter density tiq. Comparing to that of HLSinin(7'', P, ^^i 
we see that the inclusion of the 0{p'^) interaction terms and the other ir-p-uj inter- 
actions through the hWZ terms is crucial to make ni/2 larger by the factor of 1.7. 
This is quite a noticeable improvement from the previous work,'^^ where the posi- 
tion of the phase transition could not be explicitly pinned down. The higher value 
of ni/2 of HLS(7r,p, w) compared to the HLSniin(7r,p, w) result can be understood 
from the fact that the size of the single skyrmion is smaller in the former and the 
additional interactions in HLS(7r,p,w) weaken the repulsive interactions from the 
w meson. This analysis is supported by the results from HLS(7r,p). 

In addition, we find the weak dependence of E/B on the density in HLS(7r,p). 
This again emphasizes the role of the vector mesons. As we can see from Table [U 
as well as from the conclusion of Ref. Ej when only the pions are considered, the 
skyrmion carries a lot of excessive energy, while the inclusion of the p meson reduces 
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Fig. 2. The energy per baryon number is plotted as a function of L. The solid, dashed and 
dash-dotted lines are {E/ B)q( 2\, {E/ B)q, i\ and (£?/B)hWZ) respectively. 



the soliton mass a lot and almost saturates the Bogomol'nyi bound. Therefore, only 
very little interaction energy remains in the skyrmion made of pion and the p, and 
consequently it yields the (nearly) flat dashed curve for E/ B. 

In Fig. [21 we present the skyrmion energy from each term in Eq. ([T]) as a 
function of L. The plot shows that [E/ B)o^pi-^ is negative for all the values of 
density, which explicitly shows that the inclusion of the higher order terms and 
vector mesons provide attractions and lead to the weak density dependence. On 
the contrary, in the single skyrmion phase, {E / B)qij,2\ is dominant and the ratio 
{E/B)Q!^pi)/{E/ B)\y^z even increases from ~ 2 to ~ 3 as the density increases. 
However, after the half-skyrmion phase transition, {E/B)(jtp2\ starts to decrease, 
while (-B/i3)hwz starts to increase. Then, at some density the repulsive interaction 
due to the w meson becomes the dominant interaction. 



4. Remarks and perspectives 

In the present work, we discuss the skyrmion and baryonic matter properties using 
the hidden local symmetry Lagrangian up to 0{p'^) terms including the hWZ terms. 
Thanks to the holographic QCD models, the LECs can be determined by the master 
formulas and the skyrmion properties can be calculated self-consistently using two 
input parameters rup and /^. This self-consistent approach was also recently applied 
to study the bound state picture of heavy baryons.^^l 

Due to the quantitative predictive power of the present framework, we apply this 
model to investigate the skyrmion matter properties by putting the skyrmions into 
the FCC crystal. For the crystal size dependence of the skyrmion energy, we found 
that, in IILSniin(7r, p, ui) and IILS(7r, p, cj), the minimum of the skyrmion energy nmin 
appears at the same density and nmin in IILS(7r, p) appears at a very deep region. 
A noticeable result is that, in contrast to IILSinin(7r,p, w), 711/2 of IILS(7r, p, w) and 
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HLS(7r,/9) appears at density higher than the normal nuclear density uq, which is 
consistent with our expectation. 

It will be also interesting to study the role of the dilaton field in our HLS model 
and the change of hadron properties in baryonic matter. Work in this direction is 
in progress and will be reported elsewhere. 
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